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We characterize averages of YliLi \ x ~ U\ a ~ l with respect to the Selberg density, 
further contrained so that t\ G [0, x] (I = 1, . . . , q) and t\ G [x, 1] (I = q + 1, . . . , N), 
in terms of a basis of solutions of a particular Fuchsian matrix differential equation. 
By making use of the Dotsenko-Fateev integrals, the explicit form of the connec- 
tion matrix from the Frobenius type power series basis to this basis is calculated, 
thus allowing us to explicitly compute coefficients in the power series expansion of 
the averages. From these we are able to compute power series for the marginal 
distributions of the tj (j = 1, . . . , N). In the case q = and a < 1 we compute the 
explicit leading order term in the x — > asymptotic expansion, which is of interest 
to the study of an effect known as singularity dominated strong fluctuations. In 
the case q = and a £ Z + , and with the absolute values removed, the average is a 
polynomial, and we demonstrate that its zeros are highly structured. 



1 Introduction 

The Selberg density refers to the probability density function (PDF) 

1 N 
Sjv(Ai,A 2 ,A)^ ^ajl^ 

supported on t S [0, 1]", where Sm denotes the Selberg integral 

N 



SW(Ai, As, A) := j n>, Al (l-iO A2 II 

■/f0,ll n l<j<k<N 

n r(A! + i + jA)r(A 2 + i + jA)r(i + (j + i)A) 
11 r(Ai + A 2 + 2 + (7v + j-i)A)r(i + A) • 



j=0 

(see |26j for a recent review). 

Special cases of ([!]), restricted to 2A G {1,2,4} and Ai,A 2 integers or half integers, arise 
in classical random matrix theory. For example, consider an N x N random real orthogonal 
matrix. With A denoting an n\ xn 2 sub-block of elements [n\ > ri2,N > ni+n 2 ) the eigenvalue 
PDF for A T A is given by (pi with 



Ai = -(m-n a -l), A 2 = -(iV-ni-n 2 -l), 2A = 1 (3) 

(see e.g. \21\ pg.124] ). For general parameters such that ([!]) is well defined, it can be realized 
as the zero distribution of certain polynomials defined by random three term recurrences |291 
[25| [T8] , or equivalently as the eigenvalue distribution for certain matrices constructed from of 



order N (not of order TV 2 ) real numbers. For definiteness, we will interpret the variables t in 
as eigenvalues. 
The integral 

E N (n; (0,x)) :=( J f dti...dt n [ dt n+1 . . . dt N S N (X lf A 2 , A; t) (4) 

\ n J J[0,x] n J[x,i] N - n 

gives the probability that exactly n eigenvalues are contained in the interval (0, x). Let Piv(n; x) 
denote the PDF for the event that there is an eigenvalue at x, with exactly n eigenvalues to 
the left. Thus p_/v(re; x) denotes the PDF for the (n + l)-st smallest eigenvalue. It relates to the 
probabilities Q via the formula 

n , 

p N {n;x) = -^2—E N (k](0,x)). (5) 

k=0 

For its relevance to multivariate statistics, the task of computing the PDFs PN(n; x) in the 
case of the parameters ^ was undertaken some years ago by Davis [15]. His strategy consisted 
of three steps. 

1. Establish that {E]y(n; (0, x))} n= o,i,...,jv occur as part of the matrix of fundamental solu- 
tions of a certain Fuchsian matrix differential equation. 

2. To (attempt) to determine the connection matrix transforming the matrix of fundamental 
solutions implied by the Frobenius type power series solution, both about x = and x = 1, 
to the matrix of fundamental solutions involving Q. 

3. To use the differential equation to generate the Frobenius type solutions, and then from 
step 2. use the connection matrices to form the appropriate linear combination required 
to compute the series expansions of {E]y(n; (0, aj))} n =o,l,...,iV- 

Because step 2. was not successfully carried out, the strategy was never executed in its 
entirety. Here we will return to this problem, and execute all three steps as required for the 
generation of Q as a power series about x = and x = 1, and thus the corresponding power 
series of(§. 

The multidimensional integrals generalizing (|4) by the inclusion of a factor YliLi \ x ~ ti\ 2fl 
can also be characterized as solutions of a matrix Fuchsian differential equation. In particular, 
with (•) denoting an average with respect to the Selberg density 0, this allows us to calculate 
the series expansion of the averages of moments of characteristic polynomials 



ti\ 2 n- (6) 



N 

l[\:>- 

l=i 

Two applications then present themselves. The first relates to the case 2/i € Z + , when ^ is a 
polynomial. As detailed in Section 3.1 it is known that the zeros of closely related polynomials 
are highly structured. We use our ability to compute the polynomial cases of ^ to (numer- 
ically) compute the corresponding zeros for a wide range of the parameters iV, Ai, A2,A. At a 
qualitative level families of curves for the loci of zeros are identified, and for large A we provide 
an approximation to the real part of the zeros. 

The second application relates to the series expansion of ^ with fi < as x — > 1 + , inter- 
preted as an asymptotic expansion. This is a problem which has attracted previous attention. It 
turns out that then the (near) degeneracies in the spectrum of varying order are being probed, 



2 



and one has an effect known as singularity dominated strong fluctuations [6j. The latter has 
many applications in addition to the one in number theory discussed in [6]. For example, it ap- 
plies to the analysis of twinkling starlight [1] , van Hove- type singularities |5j , and the influence 
of classical periodic orbit bifurcations on quantum energy level [7] and wave function statistics 

EE]- 

Let the notation A{x) ~ B{x) mean that there exists two positive numbers C\ and C 2 
independent of x such that 

In [23J @ was expressed as a certain generalized hypergeometric function based on Jack poly- 
nomials (see e.g. |214 Ch. 13]). For the latter Yan |35j had obtained estimates of the form Q. 
Consequently it was shown in [23] (with b i->- b — 1 to account for the fact that the correction 
a 1 — y cz — I — 1 is required in the second argument of 2 .f[ i n (35)) that for \i < and x > 1, with 
\fj,\ taking on the explicit value 

2\fi\ = 2X(j - 1) + 1 + A 2 , je{l,...,N}, (8) 

one has 

N 

x~ 2flN ( II \ x - ra (! - l/ar)"°'" 1)iA log (9) 
Furthermore, for < 0, x > 1 and with such that 

2A(j - 1) + A 2 + 1< 2|/x| <2Aj + A 2 + l je {!,..., N}, (10) 



it was shown 

N 

X" 2 ^( J] |x - Xl \ 2 ») « (1 - l/^-iPH-Aa-l-a-lW (11) 
/=1 



Note that (10) implies that j on the RHS of (11) can be written 



j = int[((2|/i| - A 2 - 1)/2A) + 1], (2|/x| - A 2 - 1)/2A + 1 Z, (12) 

where int[-] denotes the integer part. 

The present study allows the results Q and (11) to be supplemented by the specification 



of the series expansion about x = 1, giving in particular the proportionalities in the leading 
asymptotic forms. 

2 The Fuchsian matrix differential equation 
2.1 Specification 

Our aim in this Section is to characterize the family of integrals 



I<f\x):= (^) J dtl ---J dtq j x dt i+^---J x dt N F(h,..., 



t N ;a;x), (13) 



where 

N 

F(t 1 ,...,t N ;a;x):=]ll} 1 (l-ti) Xa \x-ti\ a - 1 ]J \t k -t 3 \ 2X 

1=1 l<j<k<N 



3 



in terms of the solution of a Fuchsian matrix differential equation. The integrals are well defined 
for 

ReAi>-l, ReA 2 >-l, Rea > 0, Re A > 0. (14) 

In the case a = 1, 2 A = 1, which is recalled in the introduction has motivation in mathematical 
statistics, this was first accomplished by Davis |15j . 



Introduce the generalization of (13) 

dti . . . 



dt 



dt 



9+1 



dt^F(ti, . . . , i/v; a; x) 



x e p (ti — x, . . . , i/v — x) 



(15) 



where e„ denotes the elementary symmetric function 



e p (xi, . . .,x N ) :-- 



E 



3n- 



l<jl<j2<-<jp<N 



It was shown in |19| that for general parameters (15) satisfies the differential-difference equation 



(N-p)E p M(x) 



(A p x + B p )J${x) + x(x - 1)^-J${x) + D p x(x - 



where 



dx 



A p = (N- p)(Ai + A 2 + 2X(N - p - 1) + 2a) 
B p = (p- N)(\! + X(N-p-l) + a) 
A, 



P-1,9 



(16) 



p(X{N -p) + a) 
Ai + A 2 + X(2N - p - 2) + a + 1. 



(17) 



We see from (17) that the equation (|16[) is independent of g; in fact the derivation of [19] was 



carried out with q = N. However the same derivation — which relies on integration by parts 
— applies for general q = 0, 1, . . . , N, 

Let us fix q and form the column vector 



4 Q) (*) = [ J AJ(*)W,..,;v, 



(18) 



It follows immediately from (16) that satisfies a particular Fuchsian matrix differential 
equation. This in turn (as noted by Davis [15J in the special case of that work) assumes a 
simpler form if we set 



J$(z) = {x - lfH^(x), HW(x) = [H^(x)} p= o X ...,n. 
Proposition 1 We have 

where 



(19) 
(20) 



Y~ 



CTjV-l &N—2 





•• 

cri a 




4 



Y~ 



b N 
c N b N ^i 

'-. '•• 











C2 b\ 
ci b 



with 



(Ji = -B N _i = + X(i - 1) + a) 

en = -iE N _i = + A 2 + \{N + i - 2) + a + 1) 

6i = -(Ajv_i + ^N—i) +N-i = -t(A 2 + A(» - 1) + a) + JV - t 

c i = Dj V -i = (^-*)( A * + a ) ( 21 ) 



Remarks. 



1. In terms of Jg (»), the structure of the matrices on the RHS of (20) alters and the 
number of matrices increases from two to three. 



2. Suppose that instead of expanding about x = we wanted to expand about x = 1. With 
y = 1 — x and H^" (y) := Hg^(x) we see from (20) that 

dy 



y 1 - y 



(22) 



Moreover, it follows from (15) that 



p r(«) 



■1) P J 



p,N-q\ 



Ai<-»A 2 



(23) 



3. In [21] the integrals ( 15 ) with e p (t\ — x, . . . , — x) replaced by 

N 

l[(t j -l)- 1 e p (l-l/t l ,...,l-l/t p ) 

3=1 

have been considered, and a recurrence corresponding to a matrix differential system with 
a structure identical to that exhibited by (20) obtained. 

4. By examining the x — > behaviour of the integrals in (18) (see Proposition [2] below) it 
can be established that 

J{x) := det[J q a \x)] q= o,..., N ± 0. 
This was noted by Davis [15] in the special case of that work. Moreover, according to 



(19) 



J(x) := (x - l) N ( N+i y 2 H(x), 



H(x) :=det[H( Q )(x)] 9=0i .. 



.v 



while it follows from (20) that 



H'(x) 



Tr Y + + Tr Y~ )H(x). 
x 1 — x 



Thus 



H(x) = H x T * Y+ (l-x)- TrY ~ 



for some Ho independent of x, which can in fact be made explicit from the x — > 
asymptotic form. Such determinants of Selberg integrals have been considered in e.g. 
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2.2 Series solutions 



The differential equation ( 20 ) is a particular example of a matrix differential system 

d , 



dx^ 



Wo IX) 



A(x)wj(x), (j = l,...,N), (24) 



where A(x) has a simple pole at x = and is analytic for < \x\ < 1. The theory of such 
systems is well established |344 Section 24]. In particular, one has that there is a fundamental 
matrix of N + 1 linearly independent vector solutions of the form 

W(x) = U{x)x B =: U(x)e Blogx . (25) 

Here U(x) is a matrix analytic for |x| < 1 and so can be expressed as a power series while B is 
a constant matrix. If B is diagonalizable 



e 



Blogx = diag^W ..,x UN ) (26) 



where {z/j} are the eigenvalues of -B. Otherwise C _1 x slogx C has diagonal terms of the form 
x Ui while the upper triangular terms, partitioned into Jordan blocks, are polynomials in \ogx 
of degree strictly less than the row size of the corresponding Jordan block. 

In practice, unless (26) holds, the decomposition (25) is not accessed directly from (24). 
Instead one seeks solutions of the form 



oo 



Wfc(x) =x CTfc ^ Pi (logx)x', (27) 

1=0 

where the P /(logz) — which also depend on k — are polynomials in log z of degree bounded 
as for B and {<7fc}&=o,i...,jV are the eigenvalues of A\ := lim x _>o xA(x). Furthermore the {uj} in 
(26) and the must be equal modulo an integer. 
Suppose in (24) we write 

j 00 

A ^ = x!2 A i xi - 

j=0 

For a an eigenvalue of Aq, and with s = logo;, it is straightforward to show that the polynomials 
in (27) are determined by the recurrences 

p' + (<7l - ^ )P0 = 

l-l 

V , l + ((a + l)I-A )p l = ^2A l _ iPi (1 = 1,2,...). (28) 

i=0 

The simplest situation for the solution of these recurrences is when p^ = for each I = 0, 1, . . ., 
and thus there are no log terms. For this to occur it is sufficient that 

(i) Aq is diagonalizable; 

(ii) no two eigenvalues of Aq differ by an integer. 

Of these only (i) is necessary, as it possible that a pair of eigenvalues of Aq differ by an integer yet 
there is no log terms. The latter scenario occurs exactly when for some I the matrix (a + 1)1 — Aq 
is singular, the Pi for i < I are all constant vectors, and Xa=o ^i-iVi IS i n the column space of 
(<t + 1)1 — Aq. Then the second equation in (28) still permits a solution with p^ = 0. Moreover 
the solution contains an arbitrary parameter corresponding to a scalar multiple of a vector from 



6 



the null space, implying that the eigenvalues a and a + 1 of Aq give linearly independent series 
solutions. 

In the situation that conditions (i) and (ii) hold we have po is an eigenvector of Aq corre- 
sponding to the eigenvalue Ok- The vectors {pj}z=i2 ... are then given by the recurrence 



1-1 



Pi = ((a k + l)I- A^Y^Ai-m. (29) 



i=0 



Furthermore, when condition (ii) does not hold but there are no log terms, the recurrence (29) 
can be used if we first perturb the parameter Ai (say) so that (a^+l)!— Aq is no longer singular. 
After pi has been computed, we can then take the limit that Ai returns to its original value. 



2.3 Connection matrix 



Let Wfc(cc) =: [{wk(x))j]j=o,...,N be the vector of power series solutions of (20) corresponding to 
the eigenvalue of Y + . Let 



H^(x) = [(x-1)-?J$(x)] p= q_ n 



(30) 



be the vector formed from the integrals (|15j) and also satisfying (20). Form the fundamental 
matrices W(x) = [wk(x)]k=o,...,N and H^ a \x) = [Hg°' ) (x)] g= o i ... > Ar. Since the system (20) has 



dimension N + 1, and {w/ c (x)}fc =0v .. i iv and {H g (x)} 9= o,.,.,jv are bases, there exists a connection 
matrix C = \cij\i,j=v,...,N independent of x such that 



H {a \x) = W(x)C. 

We are particularly interested in the integrals in the first row of H( a \x), J ( 



which ( 31 ) gives 



(«) 



(31) 



for 



l4 a \x)] q= Q,..., N = [(w fc (x)) ]fc=o,...,JV C (32) 
With the power series of (wfc(x))o calculated according to the strategy of the previous 



section, knowledge of C will enable us to compute the power series of /< 



(a), 



Conversely, 



knowledge of the small x expansion of Iq a \x), deduced in certain restricted regions of parameter 
space, can be used to compute the entries of C. For this purpose use will be made of the Selberg 
integral ([2]), as well as a variant of the Selberg integral due to Dotsenko and Fateev [T7] 



S(p,N-p)(^l> ^2) A) := / dti . . . dtp / dtp+i . . . / dtjy 
Jo Jo Ji Ji 

N 

xji^to-ii* 3 n i*fe-*ii 2A 

«=1 l<j<k<N 



(33) 



As written, this converges for |A| 1 and —1 < Re(Ai) <0, -1 < Re(A2) and is defined 
by analytic continuation outside this region. According to |21l Prop. 4.5.1] these integrals 
satisfy the first order recurrence 



S(p,N-p)(^i> ^2, A) 



V 



sin tt(N -p+l)X sin7r(Ai + X 2 + 2 + (N + p- 2) A) 



N — p+l sin7rpA sin7r(Ai + 1 + (p — 1)A) 

x 5'(p-i,Af-p+i)(Ai, A2, A), 
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and it follows from this that ( 33 ) can be expressed in terms of the Selberg integral according to 

j 



N-p 



S( P ,N-p)(^i, A 2 , A) =<Sjv(Ai, A2, A) JJ 



j'=i 



iV-i + 1 



sinvr(A^ - j + l)Asin7r(Ai + 1 + (N - j)X) 
sin7rjAsinvr(Ai + A 2 + 2 + (2iV - j - 1)A) ' 



(34) 



Proposition 2 Let Ai,A 2 ,a,A be constrained according to (14) so that (13) is well defined 
For 

Re(Ai +a + X(q - 1)) > 



we have 

while for 
we have 



A a) {x) ~ x <T ^ g (Ai,a-l,A)S 7V _ g (Ai+a-l + 2gA,A 2 ,A) 



2/ReA < Re(Ai + a) < -2(1 - l)ReA, 



Z > 



(35) 
(36) 

(37) 



A a \x) ~ a" P ^)5^-KAi+a-l + 2/A,A 2 ,A) 

IC— S>0+ \ t — ~ 



x ^./^(AijO; - 1, A). 



(38) 



Proof Let a(x) := x(Ai + a + A(x — 1)) so that a(k) = k = 0, 1, . . . , N. The significance of 
(35) is that for this range of parameters (assumed real) 

(j = 0,...,g-l), 



a q < o-j 



u 



1,...,N). 



(39) 



To see this, we note from the form of <r(x) that (39) is equivalent to requiring a q > and thus 
Ai + a + X(q — 1) > 0. In regards to (36), for j = 1, . . . , q we change variables tj h-> xtj to obtain 



I<f\x) =X'"> ( N \ fdh... f dtqW^il-xt^l-t^- 1 H \t k -tj 

x dt q+1 ... dt N ji t^ii-uf^x-ti^- 1 n it fc -tj 

Jx Jx l=q+l q+Kj<k<N 



2\ 



2A 



q N 

x ii n i^--^i 2a - 

j=l k=q+l 



Setting x = in the integrands and comparing with ([2]) gives (36). 
The significance of (37) is that with real parameters 



j 



(j = l,...,N) j^l. 



To see this we note that for (40) to hold we require 

1--< X *<1 + - 



(40) 



(41) 



S 



where x* is the maximum of cr(x) and so 



Ai + a — A 
2A ' 



(42) 



Substituting (42) in (41) gives (37). In regards to (38), the same strategy is employed, except 



that we rescale the first q coordinates tj i— > xtj together with l — q of the next N — q coordinates, 



(this choice accounts for the combinatorial factor in (38)), and compare the resulting integrals 
with x = to both ([5} and @. 



□ 



Our strategy to use knowledge of the above asymptotic behaviour to compute the entries 
of C will require restricting attention to regions of parameter space as specified by inequalities 



(35) and (37). To then justify that our expressions are valid generally, we require knowledge of 



the anaalytic structure of the entries of C as a function of Ai in the complex plane. 



Proposition 3 Both matrices of solutions in (32) are meromorphic functions of X\ and so 
consequently is the connection matrix C . 



Proof We see from ( 28 ) that each coefficient p; in w& (x) is a rational function of the parameters 
Ai,A2,A,a with poles at the zeros of det((<Tfc + 1)1 — Aq) = 0. Furthermore, we know that 
Wfc(x) has radius of convergence R = 1 about x = 0. Hence inside this radius of convergence 



-fit 



'Wfc(x) is a meromorphic function of the parameters. 



On the other hand, consider the integrals (15). As written they are analytic functions of 



the parameters in the region (14). By suitable ordering of the integration variables 



0<ti<...<t q <x< t q+ i < . . . < t N < 1 



say, and a convention for the meaning of the power functions, the integrand is an analytic 
function of the integration variables t, and more particularly t\. 



,t q , in the region (14). 



Now let 7 be a simple contour starting at x and looping once around the origin before 
returning to x. Suppose that the integration variables t\, . . . ,t q are confined to 7 with some 
prescribed ordering. The resulting integral, Jp a q \x;j) say, is analytic in Ai. Moreover, by 
deforming 7 to loop around the origin by running along the real axis, we see that Jp a q \x; , y) is 
equal to Jp\ q (x) times a finite sum of certain phases, which are complex exponentials in Ai and 
A. Consequently Jp a q (x) itself is a meromorphic function of Ai, with poles at the zeros of the 
sum of phases. □ 

The series solutions (w q (x))o are unique only up to a scalar factor. We choose the latter so 

lim(w g (x)) (0) //( a )(x) = 1. 



that in the range ( 35 ) 



(43) 



We can now proceed to use Propositions [T] and [2] to determine the entries of C. 
Proposition 4 With C = [ck,q\k,q=o,...,N we have 

k<q (44) 
k > q. (45) 



Ck,q = 0, 

n q sin7r(A; + 1 — j)Asin7r(Ai + (k — j)X) 



sin7rjAsin7r(Ai + a + (2k — j — 1)A) 
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In particular, for k > q c^ q is periodic under Ai i— )■ X% + 1 or A i— > A + 1 and is multiplied by 
(—l) k ~ q under + Furthermore, for A an integer 



Ck,q — Ck,q 



A=0 



fc\ -p-r sin7rAi 
q) sin7r(Ai + a) 

7=1 y 



(46) 



w/iiZe for A = r/s, r and s relatively prime positive integers 



Cfc,, 



X=r/s 



Cfc mod s,ijr mods 



A=r/s 



: |fc/*J.L9/*J 



Aii— »sAi,ai— «a' 
A=0 



(47) 



Proof Suppose the parameters are such that (35) is valid. Then from (39) crfc < a q for k 
0, 1, . . . , q — 1. Taking the limit x — > in the connection formula 

JV 



(48) 



fc=0 



it follows that (44) must hold, and using too (43), that 



Suppose now the parameters are such that (37) is valid. Then (40) holds. Taking x — > in (48) 
gives that the LHS is dominated by the term k = I and thus 



a q = lim 



4 Q) (*) 



o (wi(aO)o' 



I > q. 



(49) 



But according to (43) the coefficient of x ai in (w/(x))o is the same as the coefficient of x Ul in 
r(«) 



7, (x) in the range ( 35h . Thus, using (38) in the numerator and (36) in the denominator, it 



follows from (49) that for I > q 



ci, g 



(T-qOQ S( q j_ q )(\l,a- 1,A) 

(^) Si{X,a-l,X) ' 



Now using ( 34 ) gives ( 45 ) . 



□ 



Let (wjt(x))o denote (wk( x ))o normalized so that the coefficient of x ah is unity, and thus 

oo 

(wfc(x))$ = * CTfe (l + J>,fcx fc ) (50) 



1=1 



where the coefficients pi t k are independent of x (the structure ( 50 ) holds for the generic solution 
if no log terms are present; as will be seen in Section 3.2 cases with log terms can be obtained 
as limits of the generic case) . It then follows from ( 43 ) , ( 36 ) , ( 48 ) and ( 44 ) that 

N 



4 a) ( x ) = Z)(wfc(aO)oS*(Ai, « " 1, A)5 JV _fc(A 1 + a - 1 + 2kX, A 2 , A) 

k=q ^ ' 



Cfc,, 



(51) 



Remark. 
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Let 



rO rO rl rl 

K${x) =(X -If dh... dt q / dt q+l ... dt N 

J— oo J— oo Jx Jx 

x F(ti, ...,t N ;a; x)e p (t l - x, . . . , t N - x) 



(cf. (15) and (19)), and write K^ a \x) = [^p i „(^)] P , g =o,...,A r - I n [31] the connection matrix C^ ' 1 ) 
defined so that 

H {a \x) = K {a \x)C {Q ' l) 

has been computed. 

2.4 Action of the monodromy around 

A basic question relating to the solution matrix of integrals H ^ (x) is how it transforms upon a 
circuit about the origin, x i— >■ xe 2m , or a circuit about the other singular points of the Fuchsian 



equation (19), namely x = 1 and x = oo. Consider first a circuit about the origin. Thus we 



seek a matrix Mq — referred to as the monodromy matrix — such that 



#(«)/ 



H^ a \x)M . 



(52) 



Proposition 5 Suppose the series solutions (21) have no log terms. The monodromy matrix 

(53) 



Mo in (52) is given in terms of the connection matrix C in (31) by 

M = C^DC 
where 



D = diag[e 2 ^ (Al+a - 1+(9 - 1)A) ] g= o,. 



,N- 



(54) 



Proof With W{x) specified as below (30), and pz(logx) = pi (i.e. no log terms), we see that 



W(x] 



Using this in (31) shows 



W(x)D. 



W(x)DC = H^ a \x)C- l DC. 



Comparison with (52) gives (53). 



□ 



Consider now how the solution matrix of integrals £f W (x) transforms upon a circuit about 



1. Writing y = 1 — x, H^ a \y) = H^ a \x) as is consistent with the use in (21), we seek the 



monodromy matrix M\ such that 



According to (22) 



H (a) (y) 



H {a \y) = H^{x) 



H {a \y)M 1 . 



(55) 



(56) 
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where / is the (N + 1) x (iV + 1) matrix with l's along its anti-diagonal, and O's elsewhere. 
Hence 



H {a \y) 



Ai-h-A 2 



I = (H^ a \x)M Q 



I = H {a \y) (M 



Ai«A2 



Ai«A2 



Comparing with (55) we read off that 



where the second equality follows from (52) and the third from a second application of (56). 



Mi = / M, 



AioA 2 



(57) 



For the circuit about x = 00 we set H^ a \z) : = iJW(-) and seek the monodromy matrix 



Moo such that 



H { - a \z)M OQ . 



Here no further calculation is necessary due to the fundamental relationship 

M OD M M l = I. 



(58) 



In words this says that the matrix H^ a \x) is single valued when x is traced around a contour 
starting at ^ — ie (0 < e <C 1) say, looping around the points x = 1, then 00 while not crossing 
the interval [0, 1] on the real axis and [0, —zoo) on the imaginary axis. 



We have from (53), (57), and (58) that the monodromy matrices are determined by the 



connection matrix C and the diagonal matrix D (54). The latter is periodic under any of the 



shifts Aii— >• Ai + 1, an- f a + 1 or A 1- >• A + 1, and we know from Proposition [4] that C has the 
same periodicity properties (with the qualification that each entry c^g is multiplied by (— l) k ~ q 

under a 1— > a + 1). As a consequence H^(z) 

Aj^Ai+l 

must satisfy a first order matrix difference equation 



(and similarly for ai-^a + lorA^A+1) 



H^ a \x)A 



(59) 



Aih^Ai+1 



for some (N + 1) x (N + 1) matrix A, dependent on the parameters Ai, A2, a, N and analytic 
in x. In fact in the recent work [22 the explicit form of the matrix A, by way of its Gauss 
decomposition, has been found. It turns out to be a polynomial in x. The working of |22j also 
applies to the shift a 1— > a + 1, but not A 1— > A + 1 for which the explicit form of A remains 
unknown. 

It is a basic fact that successive differentiation of the first order matrix differential equation 



( 20 ) can be used to obtain a scalar differential equation of degree N + 1 for entries of the first 
row. This scalar differential equation then has N + 1 linearly independent integral solutions 

(«)/ 



{Iq (aO}g=o,...,Ar gi ven by (13). In the case A = we have 



N 



t Xl (l-t) X2 {x-t) a - l dt 



t) X2 (t 



l dt 



N- 



\N- 



In general the scalar differential operator Cn+i of degree A^+l that has {(fi{x)) q (f2(x)) g } q =o,...,N 
as its basis, where f\ and ji are the basis of the second order differential operator C2, is called 
the iV-th symmetric power of £2 • Thus we have that the scalar differential operator correspond- 
ing to (13) with A = is the iV-th symmetric power of the scalar differential operator in the 
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case N = 1 (which is independent of A). Due to the periodicity in A, the connection matrix and 
monodromy must then correspond to this same iV-th symmetric power for all A 6 Z>q. 



Suppose A = m + 1/2, m G Z>o- For A?" odd, we see from (47) that C is the tensor product 
of its N = 1 univariate form, with the connection matrix corresponding to the (N — l)/2-th 
symmetric power as specified above but with doubled parameters. For N even, the connection 
matrix and thus the monodromy is reducible. We have that l[ a \x), 1% (x), . . . ,Ij^_^(x) span 
an invariant subspace on which the monodromy acts as the (N/2 — l)-th symmeric power of 
the univariate monodromy with double parameters, while on the quotient the monodromy acts 
as the iV/ 2-th symmetric power of the univariate monodromy. As a consequence the scalar 
differential operator will then factorize. For example, with A = 1/2, N = 2 we can check that 
the third order scalar differential operator can be written 

jjf d Ai + a X 2 + a\ 
\dx x x — 1 / 

where D is a second order operator. Here the first order operator has the solution 

x Al+Q (l - x) X2+a oc !* dt x t Xl (l - h) X2 {x - t) a ~ l f dt 2 t Xl {\ - t 2 ) X2 (t 2 - x) a ~ l \h - t 2 \ 

JO Jx 

oc/{ Q) (x), (60) 

where the first proportionality follows from the Dixon- Anderson integral [21, eq. (4.15)]. We 
remark that factorizations of differential operators play a prominent role in the recent studies 
of correlation functions in the two-dimensional Ising model [361 EZ1 [lOj [9l |8] . 



More generally, we see from (47) that for A = r/s, r and s relatively prime positive integers, 
the elements of the connection matrix C exhibit a congruence property. For N > s, n ^ s + 1 
we see that the latter is either a tensor product or block diagonal. As a consequence, in these 
cases {Ip°\x)} can be written as powers or sums of solutions of lower order equations, as with 



(60). In contrast, for generic A irrational the monodromy is irreducible and such reductions are 
not possible. 

To see this latter feature we use the fact [1] that the closure of any generic monodromy group 
contains a conjugate of the monodromy group of any Fuchsian specialization. In particular, it 
has been noted above that when A = the monodromy group is the JV-th symmetric power 
of the case N = 2. This in turn generically intersects SL 2 in a Zariski dense subgroup. Thus 
the Zariski closure of the generic monodromy group contains a copy of SL 2 in its (irreducible) 
(N+ l)-dimensional representation. In fact this subgroup of SL n + i is very nearly maximized in 
that the only permitted proper subgroups of SL 2 strictly containing it are the orthogonal and 
symplectic groups. But our monodromy matrices generically do not have eigenvalues occuring 
in complex conjugate pairs and so our monodromy group is generically irreducible. 



2.5 Computation of pN(n; x) 

With the joint eigenvalue PDF specified by ([T]) we have denoted by pn(u;x) the distribution 
of the (n + l)-st smallest eigenvalue (n = 0, 1, . . . , N — 1). According to Q, ^ and (13) this 
distribution is given in terms of {ll a \x)}k=o,..., n by 



PN(n;x) 



±^lP(x i 



S N (Xi,X 2 ,X) ^ dx 



(61) 



We can use (51 ) to compute the power series expansion of the I^\x), up to some given order, 
and thus the corresponding power series of PN(n; x). 
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0.2 



0.4 



0.6 



i' 



1.0 



Figure 1: Plot of the eigenvalue distributions p^iji] x) (n = 0, . . . , 4 reading left to right) 

1. Due to Ai and A2 being equal we have the symmetry 

b 



in the case A = 1/3, Ai = A2 
p 5 (n;x) = p 5 (A-n;x), n 



0, 1, 2. Also, ps(3; x) is given by (63) with a 



1. 



Even with an infinite number of terms the pole at x = 1 in the differential equation ( 20 ) tells 



us that generically each ijp (x) will only have a radius of convergence of unity when expanded 
about the origin. Truncating the series must necessarily reduce this value as a practical figure 



for accurately determining (61). But in keeping with (23) we have the symmetry 



p N (n; 1 - x) = p N {N - n; x) 



Ai-h-A 2 



implying that it is only necessary to calculate enough terms to ensure an accurate evaluation 
for \x\ < \. 



Another relevant point is that at in (50) is, according to (21), a quadratic function of k 



This means that a contribution of the linearly independent solution (wat(x))q, say, does not 



show itself until order x° N 



r O{N*) 



in the power series. 



We have implemented the recurrence (28) to compute the power series (50), and thus the 



power series of (51) and (61) using computer algebra software. This allows a case of (28) for 
which ((<t + 1)1 — Aq) is not invertible to be computed as the limit of a case of (28) for which 
((<t + 1)1 — Aq) is invertible. But it comes at the expense of efficiency, and restricts us to small 
values of N. Nonetheless some interesting theory can be exhibited. 
Generalizing results from [24J, it was shown in [20J that 



^ + i(ME 2/(r+1)i(r+1)w+r (/(l - x) b )) 



ME 



2(r+l),N 



(x 



(r+l)a+2r 



(1-x) 



(r+l)b+2r\ 



(62) 



Here ME^ ro (g(i)) refers to the PDF proportional to 



n^(^) n \xk-xj 

1=1 l<j<k<m 



P 



and D r+ i(-) denotes the distribution of the coordinates, assumed ordered x\ < 



< x r 



labelled by multiples of (r + 1) only. In particular, choosing r = 2 and N = 1, (62) tells us that 



p 5 (3;a:) 



1 



A=l/3 



Si(3a + 4,36 + 4,- 



-x 3a+4 (l - x) 3b+ \ 



(63) 
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Indeed we find that computation of the series of ps(3; ^)| A=1 / 3 according to the above specified 
scheme is consistent with this result. A graphical illustration of our computation is given in 
Figure 1. 



3 Properties of the average (6) 

3.1 The case \i G Z + 

The average ^ with /i € Z + , or more generally the average 

N 



v e 



(64) 



i=i 



has the special property of being a polynomial in x. Moreover, it is simply related to the inte gra l 
Ipj(x). To see this make the change of variables t\ \-¥ xt\ (1 = 1,... ,N) in the definition (13) 
to obtain 



N 



lP( x ) = x ^S N (X u a- 1, A)( J](l - xU)^ 



i=i 



where the average (•)# is with respect to the Selberg denisity Sjy(\i,a — 1, A;t). We read off 
from this that 



N 



Hi* ~ W 



r uN / i ^ 

W^(l/x) 



1=1 



5jv(Ai, A2, A) 



cn->A2+l 



cn->A2+l 



(65) 



where the second equality follows from (51) 



To gain intuition into the polynomial (64) it is useful to revise some known theory. First, a 
result of Aomoto \2\ tells us that 



N 



l\(x-t l )) = c N P i N a ^(l-2x) 



(66) 



1=1 



where cm is chosen so that the coefficient of x on the RHS has coefficient unity and 

a = (Ai + 1)/A - 1, P = (Ai + 1)/A-1 (67) 

and P N ^'^\x) denotes the classical Jacobi polynomial. This latter interpretation draws attention 



to the asymptotic N — > 00 behaviour of ( 64 ) . As shown in |3| this task is made possible by the 



existence of a duality formula, refined in [21] eq. (13.72)] to read 



N 



n(s-iO"><x(II e " l(Al 



-Aa)/2A| 



ie l \(\ 1 +\ 2 +2)/\+N-2 



1=1 



1 + e 



x 1(1 — x)e 



-Mi/ 2 - xe i6l/2 



N 



CE; 



(68) 



2/\,v 
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Figure 2: Plot of the zeros of the polynomials (64) in the cases N = 10, A = 1/3, 
Ai = A 2 = Aa, a = 3, v = 20 (left plot) and the same parameters N, a, v but with A = 3. 
The ticks on the axes are the points ±1. 



where the constant of proportionality is determined as in (66) and (-)cE a N denotes an average 
with respect to the circular /3-ensemble 



C 



n 



l<j<k<N 



The name duality formula comes from the interchange of the roles of the parameters v and iV 



on the respective sides of (68); also worth highlighting is the transformation of the coupling 



A i — y 2/A. These features characterize a wide class of duality formulas in random matrix theory 

MM- 



= l,...,vN) in ([50). 
and so (((Jn 



According to (39), for positive parameters 



We can use (65 ) to give a very efficient evaluation of (|64|). Specifically, we use (|29|) to compute 
the coefficients of pi : N {I 
ON > CTj (j = 1, ...,N- 1 
the computation of p^ jv requires / numerical matrix multiplications of N x N matrices and / 



1)1 — Aq) in (29) is always invertible. Consequently 



numerical matrix times vector operations at each step. Previously jT9j the differential-difference 



equation (|16|) has been applied to this same task, as has an order N + 1 difference equation 

The drawback of these latter two approaches relative to 



satisfied by (13) for shifts a h-> o+l [22] 



1 (x) but also I^\x) for each k 



N 



1, 



the present one is the need to compute not just / 

A consequence of our ability to efficiently compute the polynomial (64) is that we can 
investigate its zeros. As motivation, we recall that in the special case A = 1 (64) corresponds to a 



family of polynomial r-functions for the Painleve VI system [27] . The zeros of the polynomial r- 
functions for various Painleve systems have been the subject of a number of studies by Clarkson 
(see e.g. |12 [ IT U [T3]). Most interestingly, in the case of Painleve IV it was shown that the zeros 
coincided with the equilibrium points for certain vortex flows [14] . 

The free parameter A relative to the requirement A = 1 for the interpretation as a Painleve 



r-function gives certain insights into the zeros of (64) for general A > 0. The main point is that 
with 

Ai h-> Aa, A2 A6, A — > 00 
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0.05 r 



-0.05 L 



Figure 3: Plot of the zeros of the polynomials (64) in the cases iV = 10, A = 3, Ax 



Ao = Aa, 



3, v 



20 with imaginary part between —0.05 and 0.05 (zoom of right plot 



of Figure [2J super imposed with a plot of the zeros of ( 69 ) 



the Selberg density crystallizes to the zeros of the Jacobi polynomial 



P 



(a-1,6-1) 



A? 



1 - 2x) 



(69) 



(see e.g. |21| Exercises 3.6 q.5]). Consequently in this limit (64) is then proportional to (69) 
raised to the power of v. Remarkably (66), (67) tells us that in the case v = 1 this result 



persists for finite A, up to the explicit value of the parameters a, b in (69). For v > 1 and A 



finite we would expect at the very least that the u-fold degeneracy of the zeros present in the 
limit A — > oo would be broken, but nonetheless that as in the case v = 1 for there to remain 
signatures of the location of the zeros of (69) when using parameters Ai = Aa, A2 



Aft. 



Computation of the zeros shows that indeed the i^-fold degeneracy of each zero present for 
A — > 00 is broken for finite A. Instead for each of the N previously u-fold degenerate zeros 
there are now N distinctive curves in the complex plane containing the v zeros, and for large 



A these curves cut the real axis close to the zeros of (69). Because the polynomial (64) has real 



coefficients, the curves are invariant under reflection in the real axis, and this in turn tells us 
that for v odd each curve contains a zero on the real axis. The length of the curves decreases 
as A increases, in keeping with the degeneracy for A — > 00. There are also distinctive curves 
tracing groups of N zeros reading from left to right in the complex plane. At a qualitative level 



there being two classes of curves is in keeping with the duality formula ( 68 ) . The specific plots 
given in Figures [2] and [3] illustrate these general features. 



3.2 Asymptotics associated with d6 



In Q and ( |11[ ) equations implying the functional forms of the asymptotic behaviour of the 
average ^ in the limit x — > 1 + are given. By the symmetry of the Selberg density ([!]) under 
t[ 1 — y 1 — t[ and Ai <-> A2 this with x 1— > 1 — x and A2 Ai gives the functional forms of the 
asymptotic behaviour of Q in the limit x — > CP. We can use the expansion ( |5l[ ) to reclaim 
these latter behaviours, and furthermore to specify the previously unknown proportionalities. 
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The first step is to write the average ^ in terms of the integrals (13). Noting 



dt\ . . . dt N f( 



h,...,t N ) = V ( N ) [ dt 1 ...dt q f 



dt, 



valid for any / integrable on [0, 1]^ and symmetric, it follows 

1 



[x,l] N ~i 



q+1 ■ ■ 



dtNf(h, . . . , tjsr), 



N 

rife 



A 



|2ju 



5jv(Ai, A2, A) 



9=0 



a=2^i+l 



We now substitute (51) to obtain 

N 



|2/u 



1 



AT AT 



5jv(Ai, A2, A) 



g=0 fc=g 

x Sfe(Ai, a - 1, A)5Ar_A:(Ai + a - 1 + 2£;A, A 2 , X)c k , q 



(70) 



a=2/Lt+l 



By combining the inequalities (40) and (41) known from the proof of Proposition [2] with the 
leading form of (wfc(x))* as apparent from (49), the following refinement of the asymptotic 
behaviour implied by (11) is obtained. 



Proposition 6 Suppose A,Ai,A2 G M obey (14) but relax the condition on a. Rather, for a 
given I G {0, . . . , N} suppose 

(71) 



We have 



N 

rife 

i=i 



2ZA-Ai <ol< -2(Z-1)A-Ai, a = 2[j, + l. 

5,(Ai, a - 1, X)S N .i(X 1 + a - 1 + 2/A, A 2 , A) 



x| 2/x ) ~ ar 

x->0 



5jv(Ai, A2, A) 



x 1 Tj ^ Cl ' q 

' q=0 



(72) 



a=2/i+l 



At the left endpoint of (71) we have —21 A — Ai = a. According to (41) the minimum of 
o~(x) then occurs at x = I + 1/2. This means that both 07 and 07+1 are equally the smallest 
exponents and consequently 

N 

n|t 7 -x| 2 ^\ ~ lim (x a 'h(l) + x a '+ih(l + 1)) (73) 

3=1 



where x ai h(l) refers to the RHS of ( 72 ) . The limit in ( 73 ) can be computed to obtain a refinement 
of the asymptotic behaviour implied by Q. 



Proposition 7 Suppose A, Ai,A 2 G M o&ey (14) but relax the condition on a and replace it by 

a = -2l\-X 1 , I G {0, . . . , N — 1}. (74) 

We have 

1 \ 5/(Ai, a - 1, A)5jv_ z (Ai + a - 1 + 2/A, A 2 , A) 



nfe 

3=1 



x\ 2 n ~ (ki a; io L 

rr-S-OV |X| 



5 A r(Ai,A 2 ,A)r(A 1 + a + 2/A) 



7 g=0 



(75) 



a=2/x+l 
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Proof We seek the individual leading order asymptotic form of the two terms in ( 73 ) . For the 
first of these, by definition 



h (l) = (*) 1, ^ S N ^i(X 1 + a - 1 + 2/A, A 2 , A) V c lA 

Now we see from Q that the term S^v-; i n (76 1 contains factors 

r(Ai + a + (2Z + j)A) (j = 0, . . . , AT - I - 1) 



(76) 



With j = we see that this factor diverges as 1/(q + Ai + 2/A) while all other terms in (76) 
remain finite and so 



h(l) 



1 



a^-2i\-\ 1 a + Ai + 2/A 

9=0 



N\ S/(Ai, a - 1, A) gjv-KAi + a - 1 + 2/A, A 2 , A) 
/ y 5jv(Ai,A 2 ,A) T(Ai + a + 2/A) 



(77) 



In contrast to the behaviour of h(l) in the limit a —> —2/A — Ai, we see that h(l + 1), defined 
by (76) with I \— )■ / + 1 has all terms involving Selberg integrals remaining finite in this limit. 
Instead, as seen from (45), each term q j9 (q = 0, . . . ,1) diverges (but not q+i^+i = 1) and we 
have 



h(l + 1) 



N \ 5 m (Ai,a-l,A) 



a->-2U-Ai \i + 1/ Sn(Xi, A 2 , A) 



S N -i-i(Xi + a- 1 + 2(1 + 1)A, A 2 , A) 



9=0 



(78) 



The expressions (78) and (77) are closely related. First, we note from (45) that in the limit 
a — > -2/A - Ai, 

sin-7r(/ + l)Asin7r(Ai + IX) 



Q+i,. 



and 



7r(Ai + a + 2/A) sin7rA 



i i 

E c l+l,q 



(79) 



(80) 



Then we use ([2]), and the reflection formula for the gamma function to deduce that in the same 
limit 

Si+i(Xi, a — 1, A) 7r(/ + 1) sin7rA 



and 



5/(Ai, a — 1, A) sin7r(/ + l)Asui7r(Ai + IX) 
SV-J-i(Ai + a- 1 + 2(1+ 1)A, A 2 , A)r(A x + a + 2/A) 1 



Substituting (79)- (82) in (78) shows 



S' A r_ i (Ai+a-l + 2/A,A 2 ,A) 

h(l) 



N-l 



h(l + l) 

a->- 2ZA-Ai 



Now substituting ( 77 ) and ( 83 ) in ( 73 ) implies the result ( 75 ) . 



(81) 
(82) 

(83) 
□ 
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